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ABSTRACT 


The  maximum  strength  analysis  of  initially  imperfect,  axially  compressed, 
orthotropic,  sandwich  and  eccentrically  stiffened, circular  cylindrical  shells 
has  beer,  developed  through  the  use  of  Reissner's  variational  principle, 
von  Karaan-Donnell  kinematics,  and  a  deformation  theory  of  plasticity. 

For  a  given  material  and  for  the  special  cases  of  isotropic  sandwich 
cylinders  and  conventional  cylinders  with  eccentrically  located  longi¬ 
tudinal  stiffeners,  the  results  of  the  analysis  reflect  significant 
reductions  in  load-carrying  capability  in  the  range  of  "effective" 
radius-to-thickness  ratios  of  practical  interest. 


FOREWORD 


The  work  reported  herein  constitutes  e  portion  of  a  continuing  effort 
being  undertaken  at  Stanford  University  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  under  Contract  DAAJ02-68-C-0035  (Task  1F162204A17002) 
to  establish  accurate  theoretical  prediction  capability  for  the  static  and 
dynamic  behavior  of  aircraft  structural  components  utilizing  both  con¬ 
ventional  and  unconventional  materials.  Predecessor  contracts  supported 
investigations  vhicl  led,  in  part,  to  the  results  presented  in  references 
4,  15,  16,  17,  19,  a-. d  21. 
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INTRODUCTION 


The  challenge  of  designing  lightweight  structures  with  high  strength  and 

stiffness  for  aerospace  applications  has  not  yet  been  met  even  though 

significant  advances  have  been  made  in  vehicle  capability.  One  of  the 

fundamental  problems  precluding  true  structural  optimization  is  that  the 

most  common  structural  element  for  resisting  compression  loadings,  the 

thin  shell,  can  be  analyzed  for  design  purposes  only  on  a  semi-empirical 

basis.  The  history  of  the  40-year-old  enigmatic  thin-shell  instability 

1  2 

problem  through  1967  is  given  in  several  recent  survey  papers  by  Hoff  ’ 

3 

and  Stein  .  The  basic  reasons  established  for  the  discrepancy  between 
prediction  and  actual  performance  of  thin  shells  in  compression  are  the 
effects  of  prebuckling  deformations,  initial  imperfections  in  shell 
geometry,  and  boundary  conditions.  Interestingly,  except  for  noting 
the  theoretical  limitations  of  the  linear-elastic  theories  surveyed, 
neither  author  discusses  nor  assesses  the  quantitative  effects  of  inelastic 
deformations  on  the  maximum  strength  (initial  buckling  load)  of  cylindrical 
shells  in  axial  compression.  The  extreme  importance  of  inelastic  deform¬ 
ations  in  determining  the  maximum  strength  of  such  shells  has  been  brought 

4 

out  recently  by  Mayers  and  Wesenberg  .  They  determined  that,  in  the 
range  of  shell  radius-to-thickness  ratios  of  practical  interest,  signi¬ 
ficant  maximum- strength  reductions  are  obtained  relative  to  predictions 
based  upon  any  theory  restricted  to  linear-elastic  material  behavior. 

Since  thin  shells  for  resisting  compression  loads  appear  in  aerospace 
structures  in  either  stiffened  or  sandwich  rather  than  pure  monocoque 
form,  the  "effective"  radius-to-thickness  ratios  are  in  the  range  50-250. 
The  low  ratio  is  mentioned  specifically  by  Hoff*.  For  the  range  50-250, 
it  is  apparent  that  the  effects  of  initial  imperfections,  the  main 
contributor  to  theory-experiment  discrepancies  in  thin,  unstiffened 
shells  (radius-to-thickness  ratios  >  250),  must  be  minimized  and  the 
effects  of  inelastic  deformations  maximized  in  establishing  maximum 
strength.  This  maximum- strength  problem,  with  inelastic  deformations 
considered,  is  not  to  be  confused  with  the  plastic  buckling  of 
shells  treated,  for  example,  by  Lee^  and  Batterman*5.  Such  shells  are 


1 


prone  to  be  inefficient  coopered  with  shells  which,  if  devoid  of  initiel 

imperfections,  would  buckle  at  the  classical  bickling  level  of  the 

7  8 

linear  theory  arrived  at  independently  by  Timoshenko  ,  Lorenz  ,  and 

o 

Southwell  about  60  yeam  ago. 

Now,  in  view  of  the  results  of  reported  compression  tests  on  large-scale 

sandwich  and  stiffened  shells  (references  10,  11,  12)  in  the  radius-to- 

thickness -ratio  range  45-220,  maximum  calculated  stresses  Indicate  the 

presence  of  inelastic  deformations  at  maximum  load.  Therefore,  in  view 

4 

of  the  main  conclusion  of  Mayers  and  Wesenberg  and  experimental  evidence, 
the  present  study  has  been  undertaken  to  determine  the  influence  of 
inelastic  deformations  on  the  maximum  load  of  initially  imperfect,  axially 
compressed,  circular  cylindrical  shells  of  stiffened  and  sandwich  con¬ 
struction.  For  the  latter,  considerations  of  the  face-dimpling  and 
face -wrinkling  modes  of  instability  have  not  been  considered. 

The  present  study  follows  the  approach  used  in  reference  4  and  differs 

from  it  only  in  the  provision  of  appropriate  mechanisms  for  describing 

orthotropic  sandwich  and  stiffened  shells  (including  eccentricity  effects). 

That  is,  the  investigation  uses  a  modified  from  of  Reissner's  variational 
1“*  ,  , 

principle  "  in  conjunction  with  von  Karman- Donne 1 I  shell  theory  and  a 
deformation  theory  of  plasticity.  Reissner's  variational  principle  not 
only  permits  the  selection  of  the  stresses  independent  of  the  displace¬ 
ments,  but  also  facilitates  the  incorporation  of  inelastic  effects  into  the 
analysis.  Confidence  in  the  inelastic  analysis  based  on  Reissner's 

principle  has  been  established  in  reference  4  by  comparing  a  special 

14 

case  of  a  purely  elastic  solution  with  that  obtained  by  Kempner  through 
use  of  the  minimum  total  potential  energy  principle.  Similar  procedures 
were  folltwed  by  Mayers  et  al.^'^  with  respect  to  both  plates  and 
circular  cylindrical  shells.  With  confidence  previously  established, 
maximum -strength,  load-shortening  curves  have  been  obtained  for  sandwich 
shells  (isotropic  cores)  and  axially  stiffened  shells  (eccentricity 
effects  included)  of  several  different  materials.  Because  of  the 
significant  material  dependence  of  these  curves,  it  is  reasonable 
to  conclude  that  inelastic  deformations  must  be  included  in  determining 


the  maximum  strength  of  sandwich  and  stiffened  cylindrical  shall*  in 
axial  compression.  Although  only  specific  Illustrations  hava  baan  glvan 
for  ths  two  types  of  shells,  the  procedure  Is  general  and  may  be 
followed  for  any  set  of  geometric  and  material  shell  parameters  corre¬ 
sponding  to  sandwich  shells  with  orthotropic  cores  and  orthotropic 
shells  with  eccentrically  located  stiffeners. 
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GENERAL  THEORY 


STATEMENT  OF  PROBLEM  AND  BASIC  ASSUMPTIONS 

The  general  problem  studied  is  the  maximum  strength  of  orthotropic 
sandwich  and  stiffened  (eccentricity  effects  included)  circular  cylindri¬ 
cal  shells  in  axial  compression.  The  solutions  are  obtained  through  the 
use  of  a  modified  version  of  the  Reissner  principle,  the  von  Karman- 
Donnell  strain-displacement  relations,  and  a  deformation  theory  of 
plasticity.  Since  no  unloading  in  the  nonlinear,  inelastic  range  occurred 
in  the  monocoque  shell  analysis  of  reference  4,  it  is  assumed  that  the 
terms  nonlinear  elastic  and  inelastic  can  be  used  interchangeably.  For 
the  sandwich  shells,  no  provision  is  incorporated  to  allow  for  the  face  - 
wrinkling  and  face-dimpling  modes  of  instability. 

The  model  used  to  describe  stiffened  and  sandwich  cylindrical  shells  is 
shown  in  Figure  1  and  comprises  a  two-element  cros.^  section  and  core  of 
finite  transverse  shear  stiffness  in  both  the  xz-  and  yz-planes, 
respectively.  The  effectiveness  of  this  model  was  demonstrated  by 
Mayers  and  Chu  in  their  maximun-load  predictions  for  sandwich  plates. 

The  two-element  section  with  the  core  rigid  in  shear  has  been  employed 
4  15-19 

by  Mayers  et  al.  ’  in  their  extensive  maximum-strength  studies  of 

plates  and  cylindrical  shells  to  avoid  the  complexity,  due  to  inelastic 
effects,  of  integrating  a  nonlinear  stress  distribution  through  the 
thickness. 

The  strain-displacement  equations  used  in  the  present  analysis  are  those 
of  von  Karman-Donnell  for  shallow,  initially  imperfect,  circular  cylin¬ 
drical  shells.  Since  the  magnitudes  of  both  deflections  and  rotations 

4 

are  small  as  determined  by  Mayers  and  Wesenberg  ,  the  use  of  the  von 
Karman-Donnell  strain-displacement  equations  is  justified. 

BASIC  EQUATIONS 

The  von  Karman-Donnell  strain-displacement  relations,  modified  to  include 
the  effects  of  initial  imperfections,  are  well  known  (for  example,  see 
reference  4) ;  they  are  given  here  as 
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The  curvature-displacement  equations  have  been  modified  to  include  the 
effects  of  transverse  shear  (for  example,  see  reference  17)  and  are  given 
by 

J  K  ~  -w,  +  > 

X  XX  XZ  ,  X 


K  =  -w,  4  y 

V  yy  yz.y 


(2) 


.  1  ..  .  1 


kxy 


'xv  2  xz,y  2  yz,x 


The  total  strains  are  then 
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VARIATIONAL  PRINCIPLE 

Because  of  its  proven  success  in  determining  the  maximum  strength  of  plates 
and  shells  (see  references  4,  15,  and  17),  the  Reissner  principle  is 
selected  for  application  in  this  problem.  For  prescribed  surface  displace¬ 
ments  (controlled  end  shortening),  the  Reissner  functional  is  simply 
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where  F'  is  th;  stress  energy  density.  F'  is  a  function  of  the 
stresses  such  that  the  strains  are  related  to  the  stresses  by  the  rela¬ 
tionships 
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The  functional  given  by  equation  (4)  it  utilized  in  the  indirect  varia- 
tionel  approach  to  develop  Euler  equation  boundary  conditions  for  both 
orthotropic  sandwich  and  stiffened  cylindrical  shells  in  Appendixes  I 
and  II,  respectively.  For  problems  involving  linear-elastic  material 
behavior,  F1  becomes  the  complementary  energy  density.  As  shown  in 
references  4,  15,  and  18  for  a  nonlinear  elastic  material,  F'  Is  given 
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The  relationship  between  e  and  a  is  given  by  the  Ramberg- 
2o  eff  eff 

Osgood  three-parameter  representation  of  a  uniaxial  stress-strain  curve 
in  the  form 
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where  K  and  N  are  material  constants  and  depend  on  a  given  material. 
Upon  substitution  of  equations  (6)  and  (7)  into  equation  (4)  and  sub¬ 
sequent  integration,  the  Reissner  functional  becomes  _  , 
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Consideration  of  purely  elastic  core  behavior,  and  integration  over  the 
cross  section  results  in  the  expression 
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where  u>  ■  x/L,  t  ■  y/2jtR,  and  V  *nc^LR.  The  primed  end  double-primed 
quantities  refer  to  the  midsurface  and  bending  contributions,  respec¬ 
tively.  The  subscripts  t  and  b  refer  to  the  top  end  bottom  faces, 
respectively,  of  the  model  cross  section. 
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METHOD  OF  SOLUTION 


ISOTROPIC  SANDWICH. CIRCULAR  CYLINDRICAL  SHELLS 


Although  the  general  theory  has  been  established  for  sandwich  shells  with 
orthotropic  cores,  the  application  is  made  herein  to  the  special  case 
of  an  isotropic  core.  Thus,  the  following  simplifications  can  be  madel 
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Substitution  of  equations  (10)  and  (11)  into  equation  (9)  with  Poisson's 
ratio  taken  as  0.5  (incompressible  material  assumption  justified  a 
posteriori  in  RESULTS  AND  DISCUSSION)  leads  to  the  result 


11  :  a '  a ' 

-  +  -JLt'  + 

EV  JJ  {  E  x  E  y 

0G 


£"  + 
e  y 


,  G  t 
....  +  I  c  c 

/  xy  4  E  t  f 


E  x 


<4  +  v 


■  i  m  •  ph  -  w  (+)  *  •  iw 

.  m  •  <+f  -  ©  a  •  •  frri 


8 


1  K 

2  (N+l) 


(12) 


Displacement,  Stress,  and  Shear  Strain  Formulations 

The  displacements  from  the  initial  shape  are  selected  to  be 
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where  e  is  the  prescribed  unit  end  shortening.  These  displacements 

are  the  same  functions  as  those  used  in  reference  4  with  the  radial 

14 

displacement  identical  to  that  of  Kempner  .  The  rationale  underlying 
the  use  of  expressions  for  u  and  v  involving  no  free  parameters  has 
been  demonstrated  and  clearly  justified  in  references  4  and  15. 

This  significant  simplification  of  the  expressions  for  the  midsurface 
displacements  permits  the  inplane  equilibrium  equations  in  the  x-  and 
y-directions  to  be  satisfied  independently  of  the  magnitudes  of  the 
free  stress  coefficients. 

21  22 

From  the  findings  of  Mayers  and  Wrenn  and  Tennyson  a  1  Welles  ,  as 
pointed  out  in  reference  4,  an  initial  radial  imperfection  shape  to 
which  a  cylindrical  shell  Is  particularly  sensitive  is 


w  =  h 
o 


'x  r.y  2rx  \ 

*11  ^s-cas-+  52q  cos— ( 

O  X  V  o  X 


(16) 


The  midsurface  stresses  can  be  expressed  in  the  form 
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The  equivalence  of  the  coefficients  in  ~  ,  — *  ,  and  is 

E  t  t 

established  from  direct  satisfaction  of  the  inplane  equilibrium  equations 
in  view  of  the  fact  that  the  u  and  v  displacements  involve  no  free 
paramet  s . 

The  bending  stresses  are  written  as 


The  fora  of  the  shear  strains  in  the  xz-  and  yz-planes  are  taken  to  be 

compatible  with  w,  and  v,  respectively.  Both  expressions  are  then 
x  y 

written  in  the  fora 
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Elastic  Problem 

The  von  Karm^n-Donnell  strain-displacement  relations,  modified  to  include 
initial  radial  deformations  and  transverse  shear  effects  (see  Appendix  1) , 
are  substituted  into  equation  (12)  to  give 
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In  consideration  of  linear-elastic  behavior, the  bending  stresses  can  be 
related  to  the  curvatures  through  Hooke's  lav,  thus  modifying  the  Reissner 
functional  of  equation  (12)  .  This  procedure  is  followed  successfully 
and  justified  i*».  references  4,  15,  16,  and  17.  After  setting  K  *  0 
(linear  elastic  material) ,  equation  (25)  is  rewritten  as 
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This  modified  fora  of  the  Reissner  functional  leaves  only  the  displace¬ 
ments  and  eid surface  stresses  to  be  determined  by  the  variational 
procedure.  The  all-elastic,  load- shortening  curves  are  established 
by  seeking  a  stationary  value  of  the  Reissner  functional  in  equation  (26) 

with  respect  to  the  free  parameters  *  ,  A  ,  fx  ,  Ty  ,  and  o/E  (see 

”  ij  ’ij 

Appendix  III).  The  resulting  set  of  equations  is  reducible  to  four 
equations  in  the  four  -.xknowns  3/Z,  g^,  g2Q,  and  g  7.  These  equations  « 
then  used  to  construct  the  elastic,  imperfection-  ana  shear  parameter- 
dependent,  load-shortening  curves  shown  in  Figure  2  for  representative 
values  of  the  buckle  wave  parameters  4,  T],  and  6  and  the  core 
transverse  shear  stiffness  parameter  f. 

Inelastic  Problem 


Equation  (25),  consisting  of  17  free  parameters,  can  now  be  utilized  for 
the  inelastic  problem.  The  additional  parameters  over  the  four  required 
for  the  elastic  solution  appear  because  the  bending  stresses  require 
seven  additional  parameters.  Neither  these  seven  additional  free 
coefficients  a^,  a2Q,  aQ2 ,  b^,  b2Q,  bQ2 ,  and  d^  nor  the  six  mid¬ 
surface  stress  components  ,  A22*  A13’  A31’  A20’  an<*  A02  Can 
eliminated  in  terms  of  the  g^,  g^,  gg2,  and  o/E  because  of 
the  highly  nonlinear  inelastic  contribution  to  the  Reissner  functional. 
As  a  result,  the  extrema  of  equation  (25)  are  found  with  respect  to  the 
17  unknowns  by  analytic  minimization  of  the  Reissner  functional 
and  subsequent  numerical  solution  of  the  resulting  nonlinear  algebraic 
equations.  A  (R/h)  -  dependent,  load-shortening  curve  is  obtained 
and  presented  in  Figure  3  for  a  2024-T3  aluminum  sandwich  cylinder 
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with  an  isotropic  core.  The  uniaxial  stress* strain  curve  of  the  face 
sheet  material  is  shown  in  Figure  4. 

ORTHOTROPICALLY  STIFFENED  CIRCULAR  CYLINDRICAL  S.TET-T-g 

A  stiffened, circular  cylindrical  shell  is  represented  by  a  two-element 
model  with  a  core  that  is  infinitely  rigid  in  shear  (see  Figure  1).  The 
model  is  allowed  the  freedom  of  variable  inplane  and  bending  stiffnesses 
to  account  for  the  orthotropic  effects  of  stringers  ani1  rings  in  the 
longitudinal  and  circumferential  directions,  respectively.  The  general 
Reissner  functional,  represented  by  equation  (8),  with  transverse  shear 
effects  neglected  is  then  integrated  over  the  thickness  of  the  model  to 
obtain 


The  effective  stress  ir  both  top  and  bottom  faces,  modified  to  include 
variable  stiffnesses  in  the  x-  and  y-directions ,  respectively,  is 
written  in  the  form 
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Hence,  the  above  expression  ie  new  identical  to  that  used  in  references 
4,  15,  16,  18,  and  19  for  isotropic  plates  and  shells.  How,  substitu¬ 
tion  of  equation  (28)  into  equation  (27)  leads  to  the  following  result 
applicable  to  orthotropic  two-element  media  undergoing  nonlinear-elastic 
deformations: 


Correlation  of  Stiffness  Parameters 


In  reference  4,  the  thickness  t  of  the  homogeneous,  isotropic,  circular 
cylindrical  shell  is  related  to  the  distance  h  separating  the  two 
faces  of  the  two-element  model  by  equating  the  classical  buckling  stress 
of  the  homogeneous,  isotropic  cylinder  to  that  of  the  two-element  cylin¬ 
der.  For  the  orthotropically  stiffened  shell,  a  similar  procedure  is 
followed  whereby  the  classical  buckling  stresses  of  the  stiffened  shell 
(see  reference  23)  and  orthotropic,  two-element  model  are  equated  to 
establish  expressions  relating  the  bending  and  midsurface  stiffnesses  of 
the  stiffened  shell  to  those  of  the  orthotropic,  two-element  model 
(see  Appendix  IV).  The  resulting  relationships  correlating  the  stiffi.ess 
parameters  are  derived  in  Appendix  IV,  first  on  the  basis  of  zero  stiff¬ 
ener  eccentricity  effects  and  then  with  eccentricity  effects  included. 

For  the  former  case,  the  correspondence  of  stiffness  parameters  is 
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For  an  isotropic,  two-element  model,  the  bending  and  twisting  stiffnesses 
are  related  to  the  distance  h  separating  the  two  faces  and  written  in 
the  form 
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On  the  other  hand,  for  an  orthotropic,  two-element  model,  the  bending 
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stiffnesses  must  be  related  to  and  h  ,  the  distances  separating  the 

faces  normal  to  the  x-  and  y-axes,  respectively,  and  h  ,  a  fictitious 

*y 

weighted  average  of  hx  and  h^.  These  h's  are  then  related  to  the 
bendir.g  and  twisting  stiffnesses  by 
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Substitution  of  equation  (36)  into  equation  (33)  leads  finally  to  the 
set  of  relations 
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Equations  (37)  are  reducible  to  those  of  the  isotropic,  two-element  model 
used  in  reference  4  by  introducing  .';he  simplifications 

h  =  h  =  h  =  h 
x  y  xy 
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and  by  eliminating  all  quantities  subscripted  with  either  s  or  r. 

2 

The  resulting  expression  is  simply  (h/t)  *  1/3,  which  is  identical  to 

the  expression  derived,  for  example,  in  reference  4. 


The  first-order  effects  of  stiffener  eccentricity  are  established  in 
the  detailed  development  given  in  Appendix  IV.  The  fundamental  bending 
stiffnesses  appearing  in  equations  (33)  and  (83)  are  modified  approx¬ 
imately  to  give 
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where  z  and  z  are  the  distances  from  the  shell  median  surface  to 
s  r 

the  centroids  of  the  axial  and  circumferential  stiffeners,  respectively, 
and  C  is  a  nondimens ional  eccentricity  constant. 

Substitution  of  equations  (36)  into  equations  (38)  yields 
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Displacement  and  Stress  Formulations 

Equations  (13)  through  (22)  are  used  to  represent  the  displacement  and 

stress  distributions;  they  are  of  the  same  form  as  those  selected  by 
14 

Kempner  .  However,  the  nondimens ionalizing  factor  h  in  equations  (13) 
and  (16)  has  been  replaced  with  the  effective  thickness,  t,  of  the 
stiffened  circular  cylinder.  These  two  equations  are  then  rewritten  as 


,  ,  ax  fly  .  . 

+  tu  cos  -r  c°s  t  +  h 


c  as 


2flx 


(40) 

(41) 


Elastic  Problem 


The  von  K^rm/n-Donnell ,  strain-displacement  relations  presented  in 
Appendix  II  for  the  initially  imperfect,  orthotropic,  two-element  model 
are  substituted  in  equation  (32)  to  give 


0  0 


+  w,  w 


'x  o,x 


cr  ’ 

_ x 

E 


+  w,  w 

y 


°,y 


19 


1  T  * 

4  -|Z 

V,  +U.  +  W,  W.  +W.W  +W.W 

E 

x  -y  XV  X  o,y  ’y  o,x 

3"  h 
x  x 

+  -rr  -r  w, 
t  2  xx 


-  n 
y  y 


+  -r  I  W,  i  +  h  w, 

E  l  l  yyj  E  I  xy  xy 


- 1  \b  (i)2  ^  6f  -  h  fi)  (%  ♦  i-  (?*  ♦  %  $ 

-  Mi)i i)  ■ ♦  *■  (^1]  •  5  m [(^r  i^f)\  ** 

xy  t  b  J  (42) 

For  linear-elastic  material  behavior,  K  is  set  equal  to  zero,  and  the 

procedure  of  relating  the  bending  stresses  to  the  curvatures  through 

Hooke's  law  can  be  used  again,  as  in  the  case  of  the  isotropic  sandwich 

cylinder.  Equation  (42)  is  then  rewritten  in  the  form 
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Equation  (43)  now  contains  only  displacements  and  midsurface  stresses  as 
variaticnally  dependent  quantities.  A  stationary  value  of  equation  (43) 
is  then  sought  with  respect  to  the  free  parameters  A  „  ,  and  3/E 
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so|« 


(see  Appendix  V)  ;  the  resulting  set  of  nonlinear  algebraic  equations 

is  reducible  to  four  equations  in  the  unknowns  ^q2  *  an<* 

J/E.  These  four  equations  are  then  used  to  construct  the  all-elastic, 

load-end  shortening  curves  presented  in  Figure  5  for  representative 

values  of  the  buckle  wave  parameters  f.,  T),  and  1  and  the  orthotropic 

shell  stiffness  parameters  s',  s',  s'  ,  s',  s",  s",  s"  ,  and  s". 

x  y  xy  x  y  xy 


Inelastic  Problem 

In  the  inelastic  problem,  equation  (42)  must  be  used  .  As  pointed  out 

in  the  sandwich  shell  maximum- strength  analysis  of  the  present  work  and 

in  the  conventional  shell  maximum-strength  analysis  of  reference  4,  the 

resulting  number  of  free  parameters  is  17.  A  stationary  value  of 

equation  (42)  is  sought  with  respect  to  c . . ,  A. . ,  a . . ,  b. . ,  d,,,and 

ij  ij  ij  ij  *  11’ 

cr/E.  Maximum-strength,  material-dependent,  load-shortening  curves  are 
obtained  and  presented  in  Figure  6  for  2Q24-T3  aluminum  and  stainless 
steel  (1/2  hard).  The  stress-strain  curve  for  the  latter  is  shown  in 
Figure  7.  Values  of  the  buckle  wave  parameters  and  the  stiffness 


parameters  are  the  same  as  those  used  in  the  elastic  analysis. 


RESULTS  AND  DISCUSSION 


The  maximum  strength  of  initially  imperfect,  axially  compressed,  isotropic 
sandwich  and  orthotropic  stiffened  shells  (eccentricity  effects 
included)  has  been  studied  through  the  use  of  a  modified  form  of  Reissner's 
variational  principle,  the  von  Karman -Donne 11  stain-displacement  rela¬ 
tions,  and  a  deformation  theory  of  plasticity.  Maximum-strength, 
material-dependent,  load-shortening  curves  have  been  obtained  for  stain¬ 
less  steel  and  2024-T3  aluminum  cylinders.  The  results  are  compared  with 
the  load-shortening  curves  obtained  for  both  sandwich  and  stiffened 
circular  cylindrical  shells  when  the  materials  remain  linearly  elastic. 

In  conducting  the  numerical  analyses  involving  inelastic  deformations, 
Poisson's  ratios  v,  v^,  and  p^  ha,Te  each  been  assigned  the  value 
0.5.  Trial  computations,  in  which  these  Poisson's  ratios  were  varied 
between  0.3  and  0.5,  indicated  that  the  incompressible-material  assump¬ 
tion  (that  is,  Poisson's  ratios  equal  to  0.5)  relative  to  the  direction 
of  loading  led  to  negligible  changes  in  the  results  based  on  varying 

values  of  v,  v  ,  and  u  .  Poisson's  ratios  u„  and  v  have  been 
'  x  *x  y  y 

calculated  from  the  reciprocal  relations  p  *»  E  p  / E  and  v  =  D  v  /D 

y  yxx  y  yxx 

to  ensure  symmetry  in  the  orthotropic  material  constant's  matrix. 

CIRCULAR  CYLINDRICAL  SANDWICH  SHELLS 

The  effects  of  transverse  shear  and  inelastic  deformations  on  the  maxi¬ 
mum  strength  of  circular  cylindrical  sandwich  shells  with  isotropic  cores 
have  been  established  in  this  analysis.  Considerations  of  the  face- 
wrinkling  and  face-dimpling  instability  modes  have  been  neglected  since 
the  interest  herein  has  been  aimed  at  the  bending  and  buckling  in  general 
instability. 

Based  on  all-elastic  behavior,  transverse  shear  effects  markedly  reduce 
the  maximum  load  of  a  cylindrical  sandwich  shell,  as  shown  in  Figure  2. 
This  trend  is  consistent  with  the  classical  results  presented  in  reference 
24.  For  6  =  0.02  (a  reasonable  imperfection  amplitude  for  an 
"effectively"  thick  sandwich  cylinder ), the  uppermost  curve  represents 
the  all-elastic,  load- shortening  curve  for  an  isotropic  sandwich  cylinder 
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with  a  rigid  core,  r  ■  0.  The  curva  for  r  ■  i  repratanta  a  sandwich 
cylinder  with  finite  tranverse  shear  stiffness  of  the  core.  The  partic¬ 
ular  value  of  the  shear  stiffness  parameter  is  representative  of  the 
weak-core  materials  used  in  practice.  It  is  noted  that  for  r*  1,  the 
reduction  in  maximum  load  from  that  obtained  for  the  sandwich  cylinder 
with  a  rigid  core  is  6$;  thus,  for  a  larger  value  of  F  (weaker  core), 
the  reduction  in  maximum  load  would  be  substantially  greater.  A  similar 
trend  was  obtained  by  Mayers  and  Chu^  in  their  maximum  load  analysis  of 
isotropic  sandwich  plates. 

Inelastic  effects  are  incorporated  into  the  analysis  by  introducing 
actual  stress-strain  curves.  In  the  case  of  2024-T3  aluminum,  the 
uniaxial  stress-strain  curve  is  shown  in  Figure  4.  As  a  result,  maximum- 
strength,  (R/h)  -  dependent,  load-shortening  curves  are  obtained  and 
presented  in  Figure  3.  The  upper  curves  represent  all-elastic  behavior 
for  infinite  (r  =  0)  and  finite  (F  =  1)  transverse  shear  stiffnesses, 
respectively.  The  lower  curves  show  the  effects  of  inelastic  behavior 
on  two  2024-T3  aluminum  sandwich  cylinders  of  (R/h)51  200  with  T  =  0 
and  T  -  1,  respectively  As  can  be  seen,  the  maximum  load  is  reduced 
6$  by  introducing  finite  tranverse  shear  (F  =  1)  and  is  reduced 
an  additional  14$  by  including  the  effects  of  inelastic  behavior  in  the 
face  sheets.  It  can  be  noced  also  that  for  a  rigid  core  (T  =  0) ,  the 
reduction  in  load-carrying  capability,  due  to  inelastic  deformations, 
is  16$.  However,  as  mentioned  above  for  a  cere  of  finite  transverse 
shear  (T  =1) ,  the  maximum  load  was  reduced  14$  due  to  inelastic  behavior. 
It  might  be  concluded,  therefore,  that  plasticity  affects  a  rigid-core, 
isotropic  sandwich  cylinder  to  a  greater  degree  than  one  with  finite 
transverse  shear  stiffness;  this  trend  was  noted  by  Mayers  and  Chu^ 
in  their  maximum-strength  analysis  of  isotropic  sandwich  plates.  This 
phenomenon  in  both  sandwich  plates  and  cylinders  is  attributed  to  the 
fact  that  the  finite  transverse  shear  permits  stress  relief  across  each 
buckle, causing  a  lower  effective  strain  and,  consequently,  less  reduction 
in  maximum  load  due  to  inelastic  deformations.  However,  from  the  stand¬ 
point  of  inelastic  deformations  alone,  this  isotropic  sandwich  shell  of 
2024-T3  aluminum  reflects  not  only  (R/h)  -  dependence  of  the  load- 
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shortening  curve  for  a  sufficiently  low  (R/h)  ratio  (the  same  effect 
established  for  conventional  shells  in  reference  4),  but  a  significant 
reduction  in  maximum  strength. 

ECCENTRICALLY  STIFFENED.  CIRCULAR  CYLINDRICAL  SHELLS 

The  maximum  strength  of  initially  imperfect,  axially  compressed, 
eccentrically  stiffened,  circular  cylindrical  shells  has  been  obtained 
for  material  properties  corresponding  to  stainless  steel  (1/2  hard)  and 
2024-T3  aluminum. 


The  results  of  a  purely  elastic  analysis,  presented  in  Figure  5,  reflect 

load-shortening  curves  for  integral,  longitudinally  stiffened,  circular 

cylindrical  shells.  The  integrally  stiffened  shell  cross  section,  shown 

12 

in  Figure  5,  is  identical  in  geometry  to  that  used  by  Card  and  Jones 
in  their  studies  on  the  buckling  of  eccentrically  stiffened  cylinders. 

The  values  of  the  buckle  wave  parameters,  namely,  p  =  0.5  and  T|  =  0.25, 
have  been  selected  to  correspond  with  those  of  experiment  (see  reference 
12).  A  reasonable  value  of  the  imperfection  parameter,  A  =  0.5,  was 
selected;  and  the  value  of  the  eccentricity  constant,  C  =  60,  was 
adjusted  to  yield  an  accurate  representation  of  the  effects  of  eccentric¬ 
ity  on  the  maximum  load  of  both  outside  and  inside,  longitudinally 
stiffened,  cylindrical  shells  (see  reference  12) .  The  upper  and  lower 
curves  in  Figure  5  correspond  to  externally  and  internally  stiffened 
cylinders,  respectively.  Both  types  possess  the  same  integral  stiffener 
cross  section.  It  can  be  noted  that  the  externally  stiffened  shell 
carries  38fS  more  axial  load  than  its  internally  stiffened  counterpart, 
a  result  which  is  consistent  with  those  obtained  in  references  12  and 

23.  However,  this  may  not  always  be  the  case,  due  to  the  findings  of 

25 

Hutchinson  and  Amazigo  .  They  show  that  for  certain  ranges  of  the 

2 

curvature  parameter,  Z  =  — rVl-v^  ,  an  outside  stiffened  shell  is 
much  more  imperfection  sensitive  than  one  with  inside  stiffening.  How¬ 
ever,  in  the  present  analysis,  the  stiffened  shell  configuration  is  in 
the  large  Z  range  and,  therefore,  is  out  of  the  region  in  which  an 
outside  stiffened  shell  should  be  any  more  imperfection  sensitive  than 
one  with  inside  stiffening. 
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The  effects  of  inelastic  behavior  are  displayed  in  Figure  6  by  using  the 
stress-strain  curve  of  stainless  steel  (1/2  hard),  shown  in  Figure  7, 
and  perturbing  from  the  free  elastic  parameters  which  were  used  to  develop 
the  all-elastic  behavior  curves  in  Figure  5.  The  results  in  Figure  6 
show  that  for  the  externally  stiffened  shell,  the  reduction  in  maximum 
load  from  chat  obtained  in  the  purely  elastic  analysis  is  23$,  whereas 
the  reduction  in  load-carrying  capability  for  the  internally  stiffened 
shell  is  15$.  Thus,  it  might  be  concluded  that  the  effects  of  eccen¬ 
tricity  for  longitudinally  stiffened  cylinders  are  seen  to  decrease  with 

26 

the  presence  of  inelastic  behavior.  This  trend  was  noted  by  Jones 
in  his  study  of  stiffened  shells  which  buckle  initially  in  the  plastic 

range;  it  is  attributed  to  the  fact  that  the  effective  stress  is  greater 
for  outside  stiffened  cylindrical  shells  at  maximum  load  than  for 
cylinders  stiffened  on  the  inside.  This  leads  to  a  larger  effective 
strain  for  outside  stiffened  cylinders  and,  consequently,  a  greater 
reduction  in  maximum  load  due  to  inelastic  behavior.  It  should  be 
emphasized  that  despite  the  significant  reductions  in  load-carrying 
capability  due  to  inelastic  behavior,  the  integrally  stiffened  cylinders 
of  the  present  analysis  buckle  well  below  the  0.2$  offset  yield  stress 
of  the  material  (stainless  steel  (1/2  hard) ) .  Therefore,  buckling 
analyses  of  the  types  presented  in  references  5,  6,  and  26  are  not 
valid  for  application  to  the  present  maximum-strength  problem  where, 
in  the  absence  of  initial  imperfections,  elastic  buckling  would  occur. 

After  using  the  stress-strain  curve  for  2024-T3  aluminum  (see  Figure  4)  , 
it  was  determined  that  the  reductions  in  maximum  loads  of  both  outside 
and  inside  integrally  stiffened  cylindrical  shells  (with  the  geometry 
shown  in  Figure  5)  fr.im  those  obtained  in  the  purely  elastic  analysis 
were  negligible;  that  is,  essentially  elastic  behavior  governed  the 
particular  shell  analyzed.  However,  extreme  caution  should  be  exercised 
in  any  design  procedure  involving  stiffened  shell  construction  wherein 
plastic  buckling  is  precluded,  since  large  reductions  in  maximum  load 
can  occur  for  imperfect  shells  and  a  variety  of  materials  characterized 
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by  a  ■ mall  exponent  (aay,  3-5)  in  the  Ramberg-Oagood  atreas-etrain 
curve  representation.  This  behavior  is  demonstrated  not  only  in  reference 
4  but  also  by  the  previous  results  obtained  herein  for  stainless  steel 
cylinders.  In  addition,  for  larger  values  of  n  and  reduced  (R/T) 
ratios  (for  example, 75-100) ,  even  aluminum  cylinders  (N  ~  9)  will  bend 
due  to  imperfections  and  reach  a  maximum  strength  and  fail  rather  than 
buckle  plastically  from  the  undeformed  state. 

Correlation  between  experimental  and  theoretical  results  has  been 
achieved  by  comparing  the  maximum  compressive  load  of  an  outside, 
integrally  stiffened,  2024-T3  aluminum  cylindrical  shell  of  reference  12 
(L/R  =  4)  with  the  maximum  load  obtained  through  the  present  analysis. 
Since  the  stiffened  cylinder  edges  of  reference  12  are  clamped  and  those 
of  the  present  analysis  are  "effectively"  simply  supported,  the  differ¬ 
ence  between  the  classical  buckling  loads  for  clamped  and  simply  supported 
edges  has  been  added  to  the  maximum  load  of  the  present  elastic 
analysis  for  an  outside,  integrally  stiffened  cylindrical  shell.  As  a 
result,  the  maximum  load  agrees  with  the  experimental  maximum  load  of 
reference  12  to  within  2$. 

Therefore,  it  is  believed  that  sufficient  evidence  has  been  given  to 
conclude  that,  in  addition  to  boundary  conditions,  initial  imperfections, 
prebuckling  deformations,  and  eccentricity  effects,  inelastic  deformations 
should  not  be  excluded  from  initial  buckling  (maximum  strength)  analyses 
of  stiffened,  circular  cylindrical  shells  in  axial  compression.  In  scale 
structures,  the  most  efficient  eccentrically  stiffened  shell  (at  least 
for  axial  compression)  would  be  one  whose  (R/t)  ratio  precludes  plastic 
buckling,  but  is  not  so  large  for  a  given  material  that  imperfection 
sensitivity  and  inelastic  deformation  significantly  reduce  the  initial 
buckling  (maximum  load) .  A  great  deal  of  experimental  work  is  required 
on  practically  fabricated  metal  shells  of  various  materials  in  the  range 
50  <  R/t  <  250  to  discern  the  relative  effects  of  stress-strain  curve 
shape  and  initial  imperfections. 
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CONCLUDING  REMARKS 


The  maximum -strength  analysis  of  initially  imperfect,axially  compressed, 
sandwich  and  eccentrically  stiffened,  circular  cylindrical  shells  has 
been  undertaken  to  establish  the  effects  of  transverse  shear  (sandwich 
case  only),  inelastic  deformations,  and  geometry  dependence  on  the 
load-shortening  curves  for  given  materials.  Although  numerical  results 
are  given  only  for  the  special  cases  of  isotropic  sandwich  and  longi¬ 
tudinally  stiffened  shells,  the  theory  and  solution  procedure  are 
applicable  to  orthotropic  shells  in  general. 

CIRCULAR  CYLINDRICAL  SANDWICH  SHELLS 

In  the  absence  of  inelastic  behavior,  it  is  demonstrated  that  transverse 
shear  deformations  reduce  the  maximum  load-carrying  capability „ of  an 
imperfect, isotropic  sandwich  shell.  When  inelastic  deformations  are 
taken  into  account,  not  only  is  the  maximum  load  further 
reduced,  but  the  load-shortening  curves  reflect  a  significant  (R/h) - 
dependence.  This  dependence  is  not  evident  in  any  elastic  analysis. 

Future  work  should  be  directed  toward  establishing  the  effects  of 
inelastic  deformations  in  the  core  in  conjuction  with  failure  not  only 
in  general  instability  (the  present  analysis)  but  also  in  the  face¬ 
wrinkling  modes. 

ECCENTRICALLY  STIFFENED.  CIRCULAR  CYLINDRICAL  SHELLS 

In  the  absence  of  inelastic  behavior,  it  is  shown  that  outside  stiffened 
cylindrical  shells  are  capable  of  carrying  higher  compressive  leads 
than  internally  stiffened  shells.  This  result  is  consistent  with  those 
surveyed  in  references  1-3.  However,  in  the  presence  of  inelastic 
deformations,  the  present  analysis  shows  that  the  load-carrying  capability 
of  both  outside  and  inside  stiffened,  stainless  steel  (1/2  hard)  shells 
of  R/t  =  167  is  significantly  reduced.  Corresponding  calculations 
for  2024-T3  aluminum  shells  of  the  same  geometry  reveal  negligible 
reduction  in  the  maximum  load  bused  on  purely  elastic  behavior.  Thus,  cau¬ 
tion  should  be  exercised  when  dealing  with  maximum-strength  analyses  of 
shells  having  relatively  low  exponents  (for  example,  stainless  steel,  N=3) 
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in  the  Ramberg-Osgood  stress-strain  curve  formulation.  For  shells  of 
materials  possessing  higher  exponents  (for  example,  2024-T3  aluminum, 

N  *  9),  caution  should  be  exercised  as  the  (R/t)  ratio  approaches  100. 

For  such  shells,  initial-imperfection  effects  are  minimized,  but  plastic¬ 
ity  effects  are  maximized.  Until  sufficient  experimental  evidence  is 
developed  to  prove  otherwise,  it  is  believed  that  the  most  efficient 
longitudinally  stiffened  shells  for  axial  compression  loading  are  those 
which  lie  between  the  plastic  buckling  failure  mode  and  catastropic  snap 
through  failure  in  a  purely  elastic  mode.  In  addition,  it  is  shown,  on 
the  basis  of  the  stainless  steel  (1/2  hard)  results,  that  the  marked 
difference  between  the  maximum  strengths  of  outside  and  inside  stiffened 
cylinders  based  on  all-elastic  behavior  is  significantly  reduced  when 
inelastic  effects  are  taken  into  account.  As  a  result,  the  choice  of 
outside  stiffening  to  achieve  greatly  increased  elastic-analysis-predicted 
maximum  strength  relative  to  the  choice  of  inside  stiffening  should  be 
viewed  with  discretion. 

Although  limited  numerical  data  are  presented  herein,  the  overall  results 
are  consistent  with  those  established  in  reference  4  for  unstiffened 
cylinders  in  axial  compression.  That  is,  regardless  of  the  type  of 
structure  (pure  monocoque,  sandwich,  or  eccentrically  stiffened),  it  is 
obvious  that  any  initial  buckling  (maximum  strength)  analysis  of 
practically  fabricated  cylindrical  shells  must  include  the  effects  of 
initial  imperfections,  boundary  conditions,  prebuckling  deformations, 
and  inelastic  deformations.  The  material-geometry  combination  that  is 
the  most  efficient  may  be  predicted  from  application  of  the  present 
theory  and  analysis  procedure  for  sandwich  and  eccentrically  stiffened 
shells  to  a  wide  variation  in  physical  parameters.  However,  without 
substantial  experimental  effort  to  establish  confidence  in  the  analytical 
approach,  significant  weight  reduction  in  aerospace  vehicle  designs,  a 
major  portion  of  which  involve  shell  structures,  will  not  be  realized. 
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CIRCULAR  CYLINDRICAL 


Figure  2.  Purely  Elastic,  Load -Shortening  Curves  for  Isotropic 
Sandwich,  Circular  Cylindrical  Shells. 
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CIRCULAR  CYLINDRICAL 
SANDWICH  SHELL 


Figure  3.  Loau-Shortening  Curves  for  2024-T3  Aluminum,  Isotropic 
Sandwich,  Circular  Cylindrical  Shells. 
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Figure  4.  Stress-Strain  Curve  for  2024-13  Aluminum. 
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STRINGER  STIFFENED 
CYLINDRICAL  SHELL 


sR/t 


Purely  Elastic,  Load- Shortening  Curves  for  Eccentrically 
Stiffened,  Circular  Cylindrical  Shells  With  Longitudinal 
Integral  Stringers. 
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Load-Shortening  Curves  for  Stainless  Steel  (1/2  Hard) 
Eccentrically  Stiffened,  Circular  Cylindrical  Shells 
With  Longitudinal  Integral  Stringers. 


STRESS,  ksi  (TENSION  or  COMPRESSION) 


STAINLESS  STEEL 


Figure  7.  Stress-Strain  Curve  for  Stainless  Stell  (1/2  Hard). 
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APPENDIX  I 


EULER  EQUATIONS  AND  BOUNDARY  CONDITIONS  DERIVED 
FROM  A  REISSNER  FUNCTIONAL  FOR  A  SANDWICH 
CYLINDER  WITH  PRESCRIBED  END  SHORTENING 


The  Reissner  functional  for  prescribed  end  shortening  is  defined  as 

U"  =  fff  (o  e  +  a  e  +  x  y  +  i  y  +  t  y  -F')dV  (44) 

JJJ  xx  y  y  xy  xy  xz  xz  yz  yz 

V 

The  von  Karm^n-Donne 1 1  strain-displacement  relations,  modified  for  the 
two-element  cylinder  (see  Figure  1) ,  and  including  initial  radial  deform¬ 
ation  and  transverse  shear  effects,  are 

12  li 

e  =  e'  ±  e"  -  u,  +—  w,  +  w  w,  ±  —  (w,  -  y  ) 

x  XX  x  2  *x  o,x  *x  2  ’xx  'xz,x 

t,  D 


=  e'  ±  e"  =  v,  +  w,"  +  w  w,  -  —  +  —  (w,  -  y  ) 

rtjb  y  y  y  2  y  o,y  ’y  R  2  ’yy  'yz,y' 


1  2 


w  h 


xy  .  *  y*  ±  y"  *  u,  +  v,  +w,w,  +  w  w,  +  w  w. 

t,b  xy  xy  ’y  ’x  ’x  y  o,x  y  o,y  ’x 

±  h(w,  -  \  7  ~  b  7  )  (45 

’xy  2  xz,y  2  yz,x  v 

The  stresses  in  the  top  and  bottom  faces  in  terms  of  the  midsurface  and 
bending  components  are 

cr  =  cr '  +  o" 

Xt,b  x  x 

a  -  o'  ±  a"  ... , 

yfu  y  y  <4f 


T  =  T  ±  T 

xy,.  u  xy  xy 


After  integration  through  the  top  and  bottom  faces  and  the  core  of  the 
sandwich  cross  section,  equation  (44)  becomes 


L  2itR 


U"  -  t 


-//  I  e  +0g  +0g  +Og 

JJ  )xtxt  xbxb  ytyt  ybyt 
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+  t  7  +1  V  +  —  (T  y  +  T  V  )  -  (F*  +  F') 

xyt  xyt  xyfe  xy^  t ^  xz  xz  yz  yz  t  b 


£f  ci 


F 1  !  dxdy 


Substitution  of  equation  (45)  into  equation  (47)  leads  to 


(47) 


U"  *  t. 


(C i  (•■•*••>) 


1  2 
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.  ,  .12 
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(  Xyt  "  TxyJh  [“’xy  "  2  7xz,y  "  2  7yz,xj  +  t 


’yy  "  7yz,y 
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xz  xz 


+  T  7 
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-!rl<^  +  ^]  -  (Fi +  Fi)f dxdyr 


■fl^xz  '’yz J  \  -  “/l 

The  average  stresses  and  the  bending  moments  are 

o'  =  -z  (a  +  a  ) 
x  2  xt  ’S 


(48) 


o'  =  x  (a  +  o  ) 

y  2  yt  yb 


T'  =  ^  (T  +  T  ) 
xy  2  xyt  xyfe 


Mx  =  'cf  i  <v  •  V  -  -‘f h  ** 

t  b 


(49) 


H  =  -t,  ^  (o  -  o 

y  f  2  y  y 


)  =  -t,  h  o" 

>  f  y 
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(49) 


M  *  ('  -  t  )  «  t  h  t» 

xy  f  2  xyt  xy^  f  xy 

Substitution  of  equation  (49)  into  equation  (48)  gives 
L  2nR 

U"  *  2t,  ff  o'  (u,  +  y  w,2  +  w,  w  ) 

f  JJ  /  X  ’x  2  ’x  ’x  o,x' 
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12  W 
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-  F'  I 


with  cognizance  taken  of  the  fact  that 

6F  =  6a'  +  6o'  +  6^'  +  +  Sr 

x  bo'y  y  3T^y  xy  x  9My  y  xy 

The  vanishing  of  the  first  variation  of  U"  with  respect  to  cr^,  cr^, 

t'  ,  t  ,  t  ,  M  ,  M  ,  M  ,  u,  v,  w,  7  ,  and  7  requires  that 

xy  xz  yz  x  y  xy  xz  yz 
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1  1  Mx  M 
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Integration  by  parts,  as  appropriate,  leads  to 
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(52) 


For  equation  (52)  to  vanish  for  simultaneous  arbitrary  variations  in 
the  states  of  stress  and  displacement,  each  of  the  above  terms  must 
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Core -shear  cfdlibrin  equations  are 


-  X 
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lateral  equilibrium  equation  is 
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Transverse  shear  stress-strata  re lot Laos  are 
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The  aidswrface  stress-displ 


tj  integrals  are 
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The  remaining  boundary  integrals  are 


€57> 

si^lifiei, 

*y> 

1511 

(59) 
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(60b) 

(60c) 
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if  xy.y 
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Since  the  end  shortening  is  prescribed,  &  vanishes  at  x  *  0,  L.  The 
remaining  boundary  integrals,  evaluated  at  x  *  0,L,  can  be  ignored 
due  to  the  very  large  L/R  ratio  assumed  for  the  sandwich  cylindrical 
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•bell.  Iti  toriwy  cowiltlo—  «t  y  ■  0,  2x1  are  either  coatlaslty 

t 

coeditioes  or  coditlou  wfclti  preclude  ny  mvltnt  sheer  or  circaafer- 
eat  lei  load  aloof  any  geaexator  after  eabeefaeat  Integral loo  over  the 
•sell  leafth. 
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APPESDLX  II 


EULER  EQCATIOKS  A3©  SOUMiASY  COHJITICES  DERIVED  FROM 
Tiff  KE1SSXL3  FlTfCTIOBtAL  FOR  AX  OCTBOYKOPIC. 
TUO-EUDfiEXT  CTU!®tt  WITH  PRESCRIBED 
_ EXT  SFJKTEglXC _ 

The  Seissaer  function*!  for  prescribes!  cod  shortening  is 


IT"  -  fff  e  +  £  *  -  -  F5)  dV 

JJJ  XX  y  y  XV  XT 


(62) 


V 

where  the  core  of  the  tw-elneat  cylinder  has  been  assumed  rigid  in 
shear. 


The  von  Erfmln-Dunne  1 1  strain  di>l.'esent  relations  for  the  orthotropic, 
two -element  cylinder,  sonified  to  include  initial  radial  iasper feet ions, 
are 
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y  *  y1  A  y"  =  U,  +  V,  +  W,  W,  +  W  W,  +  V  V,  ±  h  W, 

J^t  b  xy  xy  *y  x  x  ’y  o,x  y  o,y  x  xy  ’xy 

where  h  ,  h  .  and  h  are  introduced  to  account  for  bending  stiffness 
x  y  xy 

orthotropy.  The  distances  h  and  h  separate  the  faces  normal  to 

x  y 

the  x-  and  y-directions ,  respectively,  and  h  is  a  fictitious  weighted 

xy 

average  of  h  and  h  .  The  orthotropic  bending  stiffnesses  are  related 

to  h  ,  h  ,  and  h  in  the  following  way: 
x  y  xy 
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Eth2 
f  *7 


4—2^) 


The  stiffnesses  D  ,  D  ,  and  D  could  oe  obtained,  for  example,  fros 
x  y  xv  r 

laboratory  tests  jn  the  particular  orthotropic  structure  of  interest. 
The  stresses  in  tens  of  nidsur face  and  bending  components  arc 


x  ~  x 


yt,b  >•  y 


r  «  Z  ’  + 

'*>'  t,b  *>'  *  xy 


After  integration  through  the  too  and  bottom  faces  of  the  two-element 
section,  equation  (62)  becomes 
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Substitution  of  equation  (63)  into  equation  (66)  gives 
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(67) 


+  >  h  -  IF*  +  F')1.  dxdJ 

yt  ^  7  f  b  1 

The  average  stresses  and  bending  aoaents  are 
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X  2  xt  *b 


*  -  13  *  O  ) 

7t  yb 


y  2  y^  y* 
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h 
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h 

M  «  t.  (T  -  T  )  *  tr  h  t" 

xy  f  2  xyt  xyfa  f  xy  xy 

Substitution  of  equation  (68)  into  equation  (67)  yields 


(68) 
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0  0 


1  2  w 

+  o'  (v,  +  r  w,  +  wf  w  -  -  ) 

y  y  2  ’y  ’y  o,y  R 
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xy  ’y  x  y  ’x  o,x  ’y  o,y  ’x 


M 

x 


w’xx  -  2^  w’yy  +  w’xy  'F'{  dxdy 


(69) 


which  is  identical  to  that  derived  in  reference  4. 


The  vanishing  of  the  first  variation  of  U"  with  respect  to  a^,  a^, 

t'  ,  M  ,  M  ,  M  ,  u,  v,  and  w  leads  to  the  Euler  equations  and 
xy  x  y  xy  » 
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boundary  conditions,  which  arc  darivcd  in  detail  In  refarence  4. 
Stress-displacement  relations: 


>F '  12 
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5  x 


X  T  2  w,x  T  W,X  o,x 


>F'  1  2 
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»F ' 

f — —  =  V,  +  U,  +  W,  W,  +  V  tf.  +  w  w, 

X  ’y  *x  y  o.x  ’y  o.y  ’x 


Moment  curvature  relations  are 
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2tc  W,xx 
x  f 
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*My  2tf  yy 


(71) 


a?'  =  .  i 

■'  ■  1  T  — —  TJ 

3M  tc  ’xy 

xy  £  J 


Midsurface  equilibrium  equations  are 


a 1 ,  +  x 1  *  0 

xx  xy , y 


+  a'  *  0 

\y»x  y.y 


(72) 


In  view  of  equation  (7.?)  ,  the  lateral  equilibrium  equation  is 

o'  (w,  +  w  )  +  o'  (w,  +  w  )  +  2  t'  (w,  -i-  w  ) 

x  XX  o . xx  y  yy  o,yy  xy  xy  o.xy' 


+  "I  +  ^  CM. 


+  M  -  2  M  )  -  0 


R  2tf  xx, x  y,yy  xy,xy' 


(73) 
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The  midsurface  stresp-d j.splacement  boundary  integrals  are 


L 

6u|  dy  =  0 
0 


(74a) 


L 

6v|  dy  =  0 
0 


(74b) 


(74c) 


T* 

xy 


6u| 


2jtR 
dx  =  0 


0 


(74d) 


The  remaining  boundary  integrals  are 
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2itR 
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•  XX  o,x  xy  y  o,y  2t  x,x 
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(75a) 
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(75b) 


2  M 


xy 
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The  boundary  integrals  at  x  =  0,  L  either  vanish  because  the  end 
shortening  is  prescribed  or  can  be  ignored  due  to  the  very  large  L/R 
ratio  assumed  for  the  orthocrop ic,  two-element  cylinder.  The  boundary 
conditions  at  y  *  0,  2?R  are  either  continuity  conditions  or  conditions 
which  preclude  any  resultant  shear  or  circumferential  load  along  ar.y 
generator  after  subsequent  integration  over  the  shell  length. 
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APPENDIX  IV 


CORRELATION  OF  THE  STIFFNESS  PARAMETERS  OF  THE  ORTHOTPOPIC 
TVO-ELEMENT  MODEL  TO  THOSE  OF  A  STIFFENED 
_ CIRCULAR  CYLINDRICAL  SHELL _ 

The  classical  buckling  stress  derived  in  reference  23  for  an  axially 
cc-apressed,  ring-  and  stringer-stiffened,  circular  cylindrical  shell  is 


where 
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Et  _ 

12(1  -  v2) 


S 


E  A 

s  s 

Etd 


(81) 


The  subscripts  s  and  r  refer  to  properties  of  the  stringers  and  rings, 
respectively.  The  signs  preceding  the  terms  linear  in  zg  and  z^ 
(distance  from  the  middle  surface  of  the  shell  to  the  centroid  of  the 
stiffener)  have  been  changed  from  positive  to  negative,  due  to  the  fact 
that  the  z  coordinate  in  reference  23  is  measured  positive  outward, 
whereas  in  the  present  study,  z  is  measured  positive  inward. 

The  classical  buckling  stress  of  the  orthotropic,  two-element  model  is 


This  equation  is  the  same  as  that  derived  by  Stein  and  Mayers  when  the 
transverse  shear  stiffness  is  infinitely  large. 

To  find  expressions  relating  the  bending  and  midsurface  stiffnesses  of 
the  orthotropic,  two-element  model  to  those  of  the  ring-  and  stringer- 
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stiffened  shell,  equations  (79)  and  (82)  are  equated.  For  z  *  z  *  0, 

8  IT 

the  coefficients  of  like  powers  of  (3  are  required  to  vanish,  thereby 
leading  to  the  relationships 


D  /2t,  n  El 

x  f  _  D  ,  s  s 

1  -  v  v  -  dD 

x  y  t 


D  /2t  -  E  I 

V  1 .  =  JL  i  .  JL1 

1  -  v  v  —  J$) 

x  y  t 


nf  V  V  l/GJ  G  J  \ 

Dxy/2tf  *  -f-  [l  -  -f  -  i  +  \  (  “Jo2  + 

i  /  E  1  E  1  \1 

1  /  s  s  ,  r  r 

2  y  y  dD  vx  AD  / 1 


2tcE  k.  Sx 

f  4 


1  =  s' 


G 

= 


:*E  k  +  (i  k„  +  |j  k 
f  3  y  2  rx  < 


where 


“r-f  [i+5+5+'s4  - v2) 

-  1  +(i  -  -2)s 

-2[l+(l+,)(i  +  i)  +  (*s) 

-  1  +(l  -  v2)  R 


k3  =  2  1  + 


)  (l  '  v2)  (‘  +  { 


k.  =  1  + 
4 


59 


Stiffener  eccentricity  effects  introduce  a  coupling  between  extension 

and  bending,  as  represented  by  the  terms  containing  zg  and  z^  in 

equations  (80),  However,  since  the  orthotropic,  two-element  model 

does  not  contain  a  corresponding  set  of  coupling  terms,  stiffener 

eccentricity  effects  are  included  by  modifying  the  bending  stiffnesses 

D  ,  D  ,  and  D  as  follows: 
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where  C  is  a  nonaimensior.al  eccentricity  constant. 
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APPENDIX  V 


METHOD  OF  SOLUTION  FOR  THE  ELASTIC,  ECCENTRICALLY 
STIFFENED,  CIRCULAR  CYLINDRICAL  SHELL 


Substitution  of  expressions  for  the  displacements  [equations  (14),  (15), 
(40) i  and  (41)]  and  stresses  [equations  (17),  (18),  and  (19)]  into  the 
Reissner  functional  given  by  equation  (43)  ,  and  subsequent  integration 
of  the  result,  yields 
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The  midsurface  stiffness  parameters  s^,  s^,  and  s^  are,  In  turn, 
related  to  their  stiffened  circular  cylindrical  shell  counterparts  in 
equations  (84)  and  (85)  of  Appendix  IV,  and  are  repeated  here  in  combined 
form  as 
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The  bending  stiffness  parameters  have  been  presented  in  equation  (39) 
of  the  text  and  are  listed  here  again  as 
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2  '  l  2  +  2(1  +  v)  (R  +  S)  +  2R  1(1  -  v2)  (1  +  v)  '  R  R*J 


The  quantities  F. .  and  H..,  defined  in  Appendix  III,  can  be  used  here 
if  h  and  7)  are  replaced  by  t  and  7]  ,  respectively. 

The  equations  resulting  from  the  vanishing  of  the  first  variation  of  equa¬ 
tion  (87)  with  respect  to  the  free  coefficients  of  the  assumed  stress  and 
displacement  functions  (equations  13-19)  are  reduced  to  four  equations 
in  the  four  unknowns  g^,  5q2>  and  °VE*  F°r  prescribed  values 

A,  |j,  7],  and  the  stiffened  shell  parameters,  a  relationship  between 
o^R  cR 

—  and  -=-  is  obtained  and  plotted  in  Figure  5  for  eccentrically 
stiffened,  circular  cylindrical  shells. 
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